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Abstract 



We evaluate analytically the cross section for Higgs production plus one jet through 
bottom quark fusion. By considering the small px limit we derive expressions for the 
resummation coefficients governing the structure of large logarithms, and compare 
these expressions with those available in the literature. 



1 Introduction 



The Standard Model (SM) [T]|2"] predicts the existence of a massive scalar particle known as 
the Higgs boson. Within this theory and its supersymmetric extensions [3], matter fields 
and gauge bosons acquire mass via the Higgs mechanism [IHE] . Although the Higgs remains 
undiscovered, experiments have placed restrictions on its mass [HE]. Supersymmetric 
theories require more than one Higgs boson. There are many more free parameters than 
in the SM, and the experimental constraints are correspondingly weaker [9"lll0]. 



The Large Hadron Collider (LHC) is expected to find the Higgs boson if it exists. To 
achieve this, various production mechanisms must be considered. The relative utility of 
each depends strongly on the Higgs' mass and couplings. While in the SM gluon fusion is 
by far the largest contribution to the total cross section, in SUSY theories with large tan f3 
bottom quark fusion can dominate, due to the enhanced bbH Yukawa coupling [TTH13] . 
Reviews can be found in Refs. [14]I15]. If we assume that the proton is composed only of 
the four lightest quarks and the gluon, the so called four flavour scheme (4FS), then the 



dominant leading order diagram for this process is that shown in Fig. 1(a) Integration 
of the phase space leads to divergences arising from the kinematical region where one 
or both bottom quarks are collinear to the initial state partons. The bottom's mass 
regulates these divergences, but they still leave traces in the form of large logarithms 
]xi(m% / m^) . Such logarithms jeopordise the convergence of the perturbative series, so 
ideally one would like to resum them. This can be achieved by introducing bottom quark 
PDFs. In this five flavour scheme (5FS) jl6|, [TT] the b quark can appear in the initial 



state, and so the leading order process is changed to that appearing in Fig. 1(b) One 
sets the b quark mass to zero in this case. Results obtained in either scheme should be 
the same, although when truncating at a finite order there will be differences, formally of 
higher order in a s . Despite this, it was found that the inclusive cross sections in the 4FS 
and the 5FS differ by roughly a factor of five when evaluated at [j>f = [J-R = ran, where 
Hf/^-r is the factorisation/renormalisation scale. This remains true also at NLO QCD 
which was calculated for the 5FS in Refs. [181119], and for the 4FS in Refs. [201121] . It was 
thus proposed in Refs. [19,22~H24] that when using the five flavour scheme the appropriate 
central scale is m#/4. Indeed, the NNLO result [25] in the 5FS seems to confirm this 
choice. 

Higgs production in association with one or more jets [26| has also received much attention. 
In the case of gluon fusion the leading order cross section is known, including the full 
top and bottom mass dependence, in both the SM [26H28] and MSSM [28]. The NLO 
corrections are known only in the heavy-top limit [29431] . As far as the MSSM is concerned, 
one expects that to a very good approximation one can simply replace the effective ggH 
coupling of the Standard Model with its MSSM value [32 34J . However, as we have stressed, 
for large tan (3 one must also include the bottom fusion contribution. That is the subject 
of this paper. 

The bottom fusion contribution to -ff+jet production has been considered for the case 
in which a final state b quark is tagged [35]. This is a useful observable because one 
can measure the bbH Yukawa coupling directly. Without b tagging, this process is a 
contribution to the total -ff+jet cross section, and must be considered alongside gluon 
fusion. For this case various distributions at NLO have been presented [36], based on 
Catani-Seymour subtraction. In this paper we give the same cross section analytically. 
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Figure 1: Diagrams for Higgs production through bottom fusion in the (a) four and (b) five 
flavour schemes. 



As well as providing a very strong check on the results of Ref. [36], our results allow us 
to analytically take the pt — > limit, and thus derive expressions for the resummation 
coefficients which govern the structure of large logarithms. This will be described in 
Section [5j 



2 Notation 

We are interested in the transverse momentum distribution of Higgs bosons arising from 
the scattering of two hadrons h\ and /12 at centre of mass energy ^fS. In particular, we 
consider in this paper only that part of the cross section proportional to the Hbb Yukawa 
coupling. The 6-quark mass is set to zero everywhere except in this coupling. For a Higgs 
boson with tranverse momentum pr and rapidity y, the cross section is 

^ ' ^dxi / dx 2 fi/ hl (xi,Hf) fj/h 2 ( x 2, Mf) lJ A , (1) 



dp 2 T dy + "~ ^V*iv~^/WV^,^, dp 2 dy > 

l 3 



where fi/fa (x%, fip) is the parton density for finding a parton i in hadron h±. We expand 
the partonic cross section appearing on the right hand side in powers of the strong coupling 
constant a s (/iR,), 



d&ij tt m\ 1 1 



dp^dy 8 V 2 s C {j 



R 



2tt 



(2) 



where V = 246 GeV is the vacuum expectation value of the Higgs field, mj, is the bottom 
quark mass, Cy is a colour averaging factor (dj = 9 for quark-quark scattering, etc.) and 
the dots stand for higher terms in the a s expansion. Although the partonic cross section 

(n) 

itself is not a function of the renormalisation scale /iR, its expansion coefficients G\j are, 
so that truncating at any finite order of perturbation theory leads to an unphysical 
dependence of the cross section. Reducing this unphysical scale dependence is one of the 
primary motivations for calculating higher order QCD corrections. 

We denote the four momenta of the incoming hadrons as Pi and P2, while those of the 
colliding partons are p\ = x\P\ and P2 = xiPi- The mass, transverse momentum and 
rapidity of the Higgs boson are written mn, Pt and y respectively. Momentum conservation 
implies 

Pl+P2 = Q+PH, (3) 

where Q represents the total momentum of the final state QCD partons, of which there 
can be either one or two. These we will label p% and p±. Our results for the coefficients 
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(n) 

G^> will be given in terms of the following partonic invariants 

s = (pi+p 2 ) 2 , 

u = (Pi-Qf, (4) 

t = {P2~Q)\ 

in terms of which momentum conservation imposes the constraint 

s + u + t = m 2 i + Q 2 . (5) 

In terms of these variables the transverse momentum of the Higgs satisfies 

2 ut - m\Q 2 
Pt = • (6) 



It is also useful to define 



Su — u — Q 2 , 

S t = t-Q 2 , (7) 

2 2,2 

m T = m u +p T , 
P T 



Q 2 +P 2 T 



At leading order (0(a s )) there are two contributing channels. We find for the correspond- 
ing coefficients = gij8(Q 2 ) with, 

g bl = AC F C A (s2+ u f R \ (8) 

g bg = 4C F C A {u2 ^ s f K \ (9) 

where Cp = 3 and Ca = 3. We will discuss how to integrate these expressions over the 
momentum fractions x\ 2 in section 



3 NLO Corrections 



Figure 2: Examples of one loop diagrams at NLO 



The 0{a 2 ) corrections, which form the coefficient function G^ , receive three different 
types of contribution: the one loop virtual corrections to the leading order processes, for 
which some sample diagrams are given in Fig. [2j the mass factorisation pieces, arising from 
the definition of the parton densities at NLO, and finally the real radiation contribution. 
Each of these pieces is divergent in four dimensions, so that in practice a regularisation 
procedure is required. We use conventional dimensional regularisation (CDR), working 
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in d = 4 — 2e dimensions, so that the divergences manifest themselves as poles in the 
parameter e. For infrared safe observables these poles cancel, and at the end of the 
calculation we may safely take the limit d — > 4. 

The one loop amplitude for bb — > Hg was given in Ref. [35], and we have independently 
checked the result. We also require bg — > Hb, which can be obtained by crossing. The 
virtual parts contain ultraviolet divergences, which we remove by renormalising a s and 
nib in the MS scheme. 

At NLO there are also additional channels to consider, beyond those which contribute 
at LO. In our case these are gg, bq, bb and qq. Here q represents one of the u, d, s or 
c quarks, and it is understood that the charge conjugated processes are also included. 
Because these channels do not contribute at leading order, none of them have a one loop 
correction. However, with the exception of qq, they do require mass factorisation to remove 
collinear poles, and so must be regularised just as in the case of bb and bg. 



3.1 Real Radiation 



-T5oT5oTT 




Figure 3: Examples of real emission diagrams at NLO. The initial state lines are bottom quarks, 
but the final state quarks can be of any flavour. 



We have evaluated the amplitudes for the Higgs plus two parton processes using FORM 
[37j . Sample diagrams are shown in Fig. [3j They are expressed in terms of the invariants 
s ij = {Pi + Pj) 2 an d Sijk = (pi +Pj +Pk) 2 ■ Note that we must retain the 0(e) pieces of the 
amplitudes, because integration over the three body phase space can generate poles in e. 
As a check one may verify that gauge invariance holds. We have also used MADGRAPH [38J 
to check the amplitudes in the e — > limit. 

We write the three body phase space factor as 

drs = VWJ VW) W? W^e) ^ dPT dy ' (10) 

and integrate analytically over the angular factor df2, given by 

[dn = —[ sin 1 - 2 ' 6 del sm~ 2t (/> dep. (11) 
J 27r Jo Jo 

It is useful to work in the Q rest frame. We can then obtain expressions [3D] for the 
invariants in terms of angles and energies. To perform the integration we first make use 
of momentum conservation and numerous partial fraction identities to ensure that each 
term in the amplitude squared contains at most two invariants that depend on the angles 
9 and 4>. These invariants are S13, S14, S23, S24, S123 and S124. For example, the relation 

1 1/1 1 

•Sl3«14 S u Vsi3 S14 
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reduces the number of angle dependent factors from two to one. Once this decomposition 
is achieved the integrals over df2 are of two types (m and n are integers): 



/ 



dn s ^sj 3 n , (is) 

dQs^s^. (14) 



A closed form result, valid for abritrary e, for integrals of the first type is given in Ref. [39J. 
The integrals of the second type are given as expansions in e in Ref. [ID]. We have used 
these results directly, supplementing them where necessary with extra 0(e) terms. 

As the integration over the momentum fractions is performed, divergences appear in 
the small Q 2 region due to terms with 1/Q 2 factors. These divergences are regulated^ by 
the Q~ 2e factor appearing Eq. ffTQj) . To expose the corresponding poles in e we use the 
distribution relation 



(Q 2 ) -1-6 -> — 5(Q 2 )A e + ( J -e -5- +0(e 2 ), (15) 



where A is the maximum value of Q 2 . The plus distributions above will appear in our 
final results. They are defined by 

[ A dQ 2 f(Q 2 ) [g(Q 2 )] = [ dQ 2 [f(Q 2 ) - /(0)] g{Q 2 ). (16) 
Jo Jo 

There are still divergences in the small pr region. These are regulated by the p7p 2e factor 
appearing in Eq. (|10p . In principle we could make use of a distribution relation similar to 
that above to expose the divergences as poles in e. We would then need to add the two 
loop corrections to the process bb — > H, as well as some extra terms arising from mass 
factorisation. The result would be a NNLO result for the differential cross section of the 
process bb — > H. We do not take this extra step, which is difficult to achieve in practice, 
so that in our numerical results we must avoid the small pt region. Our results therefore 
constitute a NLO result for the process bb — > H + jet. We will discuss the behaviour of the 
cross section in the pt — > limit in more detail in Section 



4 Results 



Our result for the NLO bb coefficient function takes the form 



Gf = D bl 5(Q 2 ) + E bl f±) + F bl ^ In ^) + H bl . (17) 



We define 

1 



« = ^{rn 2 K + s-Q 2 ), (18) 

A 2 = K^-m^s, (19) 

* = =±£ (20) 

K — A 



1 Strictly speaking, we mean 'made integrable'. 
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and introduce the following convenient logarithm abbreviations. 



If 





s 










-t 


In 




L u 


i ~ u 


if 


= ln 






m H 






m\ 


In 


l4 

2 ' 


Ir 


= lnii 




= ln 


A 




m H 













(21) 



We give here the D bb , E bb and F bb coefficients. The expressions for H bb are very large, so 
we do not reproduce them here. Results for all the above coefficients, for all channels, are 
available from the author on request. 



D. 



66 



9bb C A 



L A — 2LaL s + 2L A L t + 2LaL u 



11L A , lll R 



i. 



+ 



6 



Li + 2L s L t 



+ 2L S L U — L 2 — 2L t L u — L 2 U — 2Li2 ( 1 



mi 



+ 



67 
18 



+ 9bb C F 



AL 2 A + AL A L S - AL A L t - AL A L U - \L A l F + 3l R + 2L 2 S 



AL s L t - 4L S L U + L 2 t + 2L t L u + 2L t l F + L 2 U + 2LJ F - 3l F 



+ 4Li 2 1 



mi 



+ In 2 1 



mi 



+ 2Li 2 



+ In 2 1 



mi 



u 

2 

mi, 



+ 



-n f T R g bb 



L A -Ir- - 



+ 2Li 2 

^ 2 

2 

3 

1 1 



mi 



4 



+ 4m&( - + -) C F C A (C A -C F ), 



66 



+ ^66 



ln(l _ v ) + l n j» + In ^ + 2 ln ^ - L„ 



mi 



2 In ^ - 2 ln ^ - 4Z F - 4 ln ^-±^1 + 2 L £ 



+ 4C f Ca C f 



u t 
C A \ lncc 



4 4 

2_» + __h 



m 2 ^ 



6 ; 



+ ^T R n F g bb 



2 , s 

m H + - 



it 

s + — + u 



(22) 



(23) 



(24) 



Recall that S u , St and v were defined in Eq. (J7|), while g bb was given in Eq. (jHJ). Li 2 denotes 
the dilogarithm function. For QCD the colour factors take the values C F = I, C A = 3 and 



l 

2' 



To integrate the distributions above over the PDFs, it is useful to arrange for Q 2 to be 
one of the integration variables. This is easily achieved [H], with the result that we can 
replace 



dxi / dx 2 6{Q Z 



dxi 



x\ - xu Jo 



dxs 



X2 ~ X T Jo 



dQ 2 , (25) 



6 



where 




y 



x± = 



rriT+PT ± y 



7^ e 



A 2 



Ai 



xi(l - xt) - xu + — , 

X + (X 2 ~ X T ) ~ X V X2 + — . 



The expression for A\ or A2, as appropriate, should be substituted in Eq. (|2"2"j) in place 
of A. In this way one can obtain numerical results for the hadronic cross section. Phe- 
nomenological analyses have already been presented in Ref. [36J, so we do not repeat it 
here. 

4.1 Checks 

We have performed a number of checks on our results. Firstly, the dependence of the cross 

section on //p and /^r can be predicted due to the requirement that physical observables 

(n) 

must be independent of these scales. Our expressions for the perturbative coefficients G\j 
satisfy these constraints. Secondly, as we will discuss in the next section, the small pt 
behaviour can be compared to known resummed formulae. 

The strongest check is a comparison with a Monte Carlo numerical code [36] based on the 
Catani Seymour [12] subtraction formalism. We find excellent agreement for all channels. 
In the case of the qq channel, one can also compare to the known total cross section [25] 
by numerically integrating Eq.([T]) over pt and y. 

5 Small pt limit 

For small pt the convergence of fixed order perturbation theory is spoiled at any finite 
order by terms of the form 



To obtain reliable physical predictions for observables in this region one must resum these 
enhanced terms. The method to perform this resummation is known. The result [43] is 
expressed as an integral over the impact parameter b, 




n = 0,1,2,3. 



(26) 




(27) 
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where the bottom quark mass mi,, implicit in the prefactor do, is evaluated at the scale 
fj,R = mn- The Sudakov form factor W(b) contains the large logarithms, and is defined as 



W(b) = (C H (a s (b /b)) ® fi) b/b ) {C- bj (a s (b /b)) ® /,) (x° 2 ; b/b ) 



x exp 



A{a s (q))\n^ + B(a s {q)) 



(28) 



where <8> indicates convolution, partons i,j are implicitly summed over and 60 = 2e 7B , 
with Euler's constant. The PDFs fi and fj are evaluated at the scale b/bo, and 



x l,2 



(29) 



The resummation coefficients A, B and Cy can be expanded perturbatively, 



A(a s ) 
B{a s ) 



n=l 
oo 



n=l 



\2ir 



n=l 



C l3 {a s ) = 8 ij 5(l-z) + J2(¥:) n C,, 



2tt 



(30) 
(31) 
(32) 



The coefficient A^ 1 ) controls the leading logarithmic (LL) terms, while A^ 2 \ B^ and 
C^j give the next to leading logarithmic (NLL) terms, etc. They can be evaluated by 
performing a fixed order calculation and comparing to the resummed expression. 



5.1 Fixed Order Expansion of the Resummed Formula 



One cannot naively expand Eq. ([27)) in powers of a s . Instead, we first integrate by parts 
(we can ignore the surface term) to obtain 



da 



1 



dp^dy 



2S 



Pt 



dW(b) 



bdbJ^b) 
at) 



(33) 



With the p\ pole now manifest, we can expand. We use the DGLAP equation to evolve the 
PDFs to an arbitrary scale /Up. We also evolve the QCD coupling a s and the bottom quark 
mass rrib from their values at the scale mn or q to an arbitrary scale /Ur. The expanded 
resummed cross section is, using the notation of Ref. [3Tj, 



da 
dp^dy 



5 p 2 T 



2 2m-l 

E E 

m=l n=0 



V27T 



xC„ In 



+ 0{al) 



(34) 
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The coefficients m C n are related to the resummation coefficients as follows, 



1C1 
iCo 

2C3 
2C2 



iCx 



2 Co 



^ (1) / fe (x?)/ 5 (^) + {p bl ® fi)(4)fi(4) + fb(4)(Pii ® /i)(*2). 



l r 
2 
3 
2 



A (i) 
^(i) 



M4)h(4), 



ft _ 2BW - AW In i^- (P u M4)fb(4) + AW(CW ® 



-(P&i 8> ® /■,•)(*§) - (P 6i Pij fj){4)fb^ 



(i) 



A(*lU(*2) 



-A( 2 )-^ (1) -^ (1) ln i % 



In 



to?, 



+ 



mf 



(Ai®/<)(aj?)/s(S2) 



(i) 



+ 



C.3 



A W " + ^(2) + ^/ B (D ln /% 



2 2' 

- /3 )(c fe ( i 1) ® M4)fb(4) + (P b f ® /<)(s?)A(s8) 

o 6, £2}. 

where Cn is the Riemann ("-function ((3 = 1.202. . . ) and /3o = (11CU — 4Tr?t,j)/6. The 
corresponding expansions for Drell Yan production [33] and Higgs production through 
gluon fusion [31] have been presented before. Our case is slightly different to each of these 
due to the /^r dependence of the Hbb Yukawa coupling. This is reflected in the modified 
beta coefficient f3' Q = /3q + 3Cf- The two loop splitting function can be extracted 
from the results of 145,36]. 



5.2 Extracting the Resummation Coefficients 

We have checked analytically that in the limit of small Higgs transverse momentum our re- 
sults reproduce the resummed result, when the latter is expanded to the appropriate order 
in a s . Taking the limit analytically requires great care, because as well as explicitly sin- 
gular terms appearing in our results, some logarithms of px appear only upon integration 
over the momentum fractions x\p. 

By comparing with Eq. (|34p we can derive the values of the resummation coefficients A, 
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Figure 4: Diagram controlling the process-dependent part of the NNLL resummation coefficient 
B^ 2 \ as described in Ref. @S]. 



B and Cjj. For the universal coefficients we find the expected values, 

A&> = 2C F , (35) 
^ = 2C F (^C A -fn f T R -^C A y (36) 
tfW = -3C F . (37) 
For the process-specific coefficients we have 

3 o \ „ / 61 11 , \ /5 4 



# (2) = C| (-jj + vr 2 - 12Ca) + C F C A (-§ + y^ 2 + 6Cs) + C F T R n f 



vr 2 

3 9 



i-.r+(y-l),)(l-.r) 



(38) 



= 2T fl x(l-x). 



The expressions for C^- match those given in Ref. |47| . Our result for 5( 2 ) is the first 



direct calculation of this quantity. It has been shown |48j that B^ can be split into 
universal and process dependent parts, and that furthermore, the process dependent part 
is directly related to the finite part A of the one loop correction to the leading order 
process, which in our case is bb — > H. For quark initiated processes, the relationship is 
expressed as 

B® = -2 7 < 2 ) + A, (~CW 2 + A\ , (39) 
where 7 ( 2 ) is the coefficient of 5(1 — z) in the two loop splitting function P^\z), given by 



7 (2) = Cl ( 



qq 

+ ^ {I + Ts* 2 ~ 36 ) " C ™ T « (l + V I ' < 40 » 



It is straightforward to evaluate the one loop correction to bb — >■ H (the single contributing 
diagram is shown in Fig. HD, from which we find 

A = C F (-2 + ^A. (41) 
Substituting this into Eq. (|39p yields the same expression for B^ 

as we have derived from 

our analytic form of the cross section. 
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6 Conclusion 



We have described the analytic calculation of the cross section bb — > -ff+jet. The partonic 
cross section is a distribution in Q 2 , the invariant mass of the final state QCD partons. The 
results agree numerically with an implementation based on Catani-Seymour subtraction 

By taking the limit of small Higgs transverse momentum, we have evaluated the resum- 
mation coefficients that govern the structure of large logarithms, including the NNLL 
coefficient B^ 2 \ This is the first direct calculation of this quantity for this process. It 
agrees with the general expression |48j relating to a one loop amplitude. 

As well as being of phenomenological interest in their own right (numerical analyses have 
already been presented [36])) our results can form part of a differential NNLO calculation, 
perhaps along the lines of Ref. |49j . 
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